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ABSTRACT 

We assess the accuracy with which D,m and J^a can be measured by combining various types of 
upcoming experiments. Useful expressions for the Fisher information matrix are derived for classical 
cosmological tests involving luminosity {e.g., SN la), angular size, age and number counts. These geo- 
metric probes are found to be quite complementary both to each other and to inferences from cluster 
abundance and the cosmic microwave background (CMB). For instance, a joint analysis of SN la and 
CMB reduces the error bars on IIa by about an order of magnitude compared to a separate analysis of 
either data set. 

Subject headings: galaxies: statistics — supernovae: general — large-scale structure of universe — 
CMB 
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1. INTRODUCTION 

It may be possible to measure cosmological parameters 
with great accuracy using upcoming cosmic microwave 
background (CMB) experiments (Jungman et al. 1996; 
Bond et al. 1997; Zaldarriaga et al. 1997), galaxy surveys 
(Tegmark 1997; Goldberg & Strauss 1998; Hu et al. 1998) 
and supernova la (SN la) searches (Goobar & Pcrlmutter 
1995; Perlmutter et al. 1998; Garnavich et al. 1998). How- 
ever, no single type of measurement alone can constrain 
all parameters, as it will inevitably suffer from so-called 
degeneracies in which particular combinations of changes 
in parameters leave the result essentially unaffected (Bond 
et al. 1994, 1997; Zaldarriaga et al. 1997; Metcalf & Silk 
1998; Huey et al. 1998). Fortunately, different types of cos- 
mological measurements are often highly complementary, 
breaking each other's degeneracies and combining to give 
much more accurate measurements than any one could 
give alone. For example, CMB measurements are highly 
complementary to both galaxy surveys (Tegmark et al. 
1997; Hu et al. 1998; Gawiser & Silk 1998; Webster et al. 
1998; Eisenstein et al. 1998) and SN la (Zaldarriaga et al. 
1997; Tegmark 1997; White 1997). 

The topic of this Letter is probes of the acceleration of 
the Universe, given by the density parameters fim for mat- 
ter and r^A for vacuum density (cosmological constant). 
Most of the cosmological tests that we discuss are well- 
known. Our focus is on their degeneracy structure, i.e., 
on which ones are complementary and which ones act as 
independent cross-checks of one another. We address this 
by computing the Fisher information matrix F for each of 
the tests. This has the advantage of explicitly showing how 
the accuracy and degeneracy depends on the survey de- 
tails. It also allows a unified treatment of all tests, since if 
independent experiments are analyzed jointly, their Fisher 
information matrices simply add. 
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2. CALCULATION OF THE FISHER MATRICES 

All data sets discussed below consist of a vector x of 
measured numbers xi,...,xn whose probability distribu- 
tion /(x;0) depends on a vector of cosmological param- 
eters that we wish to estimate. In our case, 9i = f2m 
and 02 — ^A- The Fisher information matrix for a data 
set (see Tegmark et al. 1997 for a comprehensive review), 
defined as 

de.dOj / ' ^ ' 

quantifies its information content about these parameters. 
Its inverse F^^ gives the best attainable covariance ma- 
trix for the measurement errors on these parameters, il- 
lustrated by the error ellipses in Fig. 2. We will now spec- 
ify probability distributions / for the various cosmological 
tests and compute the corresponding Fisher matrices. 

2.1. Luminosity, size, age and clustering 

The cosmological tests based on luminosity, angular size, 
age and clustering (see e.g. Weinberg 1972, hereafter W72; 
Peebles 1993), can all be described as noisy measurements 
of some quantities Xn at redshifts z„, n = 1,...,N. We 
model them as 



(2) 



where a and b are constants independent of Sl^ and 17 a, 
the function d incorporates the effects of cosmology and 
£n is a random term with zero mean ((£«) = 0) including 
all sources of measurement error. 

For luminosity tests like SN la, Xn is the observed mag- 
nitude of the n^^ object and d is the luminosity distance 
(W72): 



dlum = (1 + z) ^^'^^^ , ri{z;nra,ftA) = 



dz' 



Eiz'Y 

E{z) = [(1 + zf[l + n„,z) - z{2 + z)nA] ^'^ , 



(3) 



(4) 



where k = a/ 1 1 — ri„i — ^a \ . We recognize E = H{z)/Ha 
as the relative expansion rate at an earlier time and 1/Hqk 
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as (the magnitude of) the current radius of curvature of the 
Universe. From the definition of magnitudes, a = 5/ In 10. 
The errors e„ include errors in extinction correction and 
intrinsic scatter in the "standard candle" luminosity. 

For tests involving the observed angular sizes 6'„ of ob- 
jects at redshifts zi, zn, we define x„ = ln0„, a = —1 
and take d to be the the angular size distance (W72): 
dang = dium/{^ + z)^ ■ For such tests (see e.g. Daly 1998; 
Pen 1997), £„ includes scatter in the "standard yardstick" 
size. 

For tests involving estimates i„ of the age of the Uni- 
verse at redshifts z„, we define Xn = InTJo^n- Setting 
a = 1, this gives (W72) 



dz' 



{l + z')E{z') 



(5) 



For tests involving the observed growth G„ in the am- 
plitude of linear density fluctuations since redshift z„, we 
choose a:„ = — In G„ , a = 1 and take d to be the linear 
growth factor (W72): 



D{z) cx E{z) 



(1 + ^) 
E{zY 



dz. 



(6) 



Assuming that the errors have a Gaussian distribu- 
tion, the Fisher matrix is given by (Tegmark et al. 1997) 



= ^tr [C- 



^C„C-iC„]+/x,*C-V,,, (7) 



where /x = (x) is the mean [/i„ = hid^Zn)] and C = 
(xx*) — /x/x* is the covariance matrix of x. Commas de- 
note derivatives, so /x,i = dfi/dOi. For simplicity, we will 
assume that 

^7nn — ^mn^n') (^) 

i.e., that all the magnitude errors e„ are uncorrelated. Our 
treatment below is readily generalized to non-diagonal er- 
ror models C, more appropriate for describing systematics. 
Since C,i = 0, all the information about 17™ and JIa comes 
from the second term in equation (|^), giving 
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2.2. A supernova example 



To bring out the physics, let us evaluate this explicitly 
for the SN la example — the other cases are analogous. 
SN la have had their accuracy assessed previously, first by 
Goobar & Perlmutter (1995) and subsequently by making 
X^-fits to real data (Perlmutter et al. 1998; Garnavich et 
al. 1998; White 1998); however, this is the first treatment 
involving their Fisher matrix. 

In this illustration, we take all magnitude errors to be 
equal, ct^ — Am. It is instructive to rewrite equation (|9|) 
as 

N r°° 

g{z)wi{z)wj{z)dz, (10) 



where 



Wi{z) 



(Am) 



KS'[Kr]{z)] 
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(11) 



d7] 1 



drj 



z'(l + z')^ 



9f]A 2 



^(^0 
z'{2 + z 



■dz' , 



dz\ 



(12) 
(13) 



E{z')- 

and the SN la redshift distribution is given by g[z) - 
W ^(-^ ~ -^n)- The expression in braces approaches 
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FIG. 1 — The SN la weight functions uia (positive) and Wm 
(negative) are plotted for standard CDM, two open (Qa = 0) mod- 
els and two flat {Q/^ = 1 — Qm) models. The Fisher matrix element 
Fij is computed by simply integrating the product of the curves Wi 
and Wj and a redshift distribution / such as the shaded one. 
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FIG. 2 — 68% confidence regions are shown for the upcoming 
CMB experiments and hypothetical SN la data sets specified in Ta- 
ble 1. The assumed fiducial model is COBE-normalized ACDM with 
f2m = 0.35, Ha = 0.65, Qt = 0.05, and h = 0.65. Combining the 
CMB and SN la data shrinks the error region to the overlap of the 
two corresponding ellipses: for instance, a joint analysis of the op- 
timistic SN la case with polarized Planck data gives the tiny black 
ellipse in the center. 
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r]^^dr]/d9i — 77^/6 as k — > 0. The contribution to F from 
each redshift can thus be spht into two factors, one re- 
flecting the quahty of the data set {N g[z] / Arn^) and the 
other incorporating the effects of cosmology (the weight 
functions Wi). The functions Wi are plotted in Figure 1 for 
a variety of cosmological models. 

If all the observed supernovae were at the same redshift 
z, then the resulting 2x2 Fisher matrix cx Wi{z)wj{z) 
would have rank 1, i.e., be singular. The vanishing eigen- 
value would correspond to the eigenvector {wq,—wa). 
Physically, this is because there is more than one way 
of fitting a single measured quantity diumiz) by varying 
two parameters ($7™ and I^a). The corresponding ellipse 
in Figure 2 would be infinitely long, with slope —wn/wA, 
the ratio of the magnitudes of the ilm and curves in 
Figure 1 at that redshift. The SN la ellipses plotted in 
Figure 2 correspond to a range of redshifts, with / being 
a Gaussian of mean z and standard deviation Az given by 
Table 1 . This breaks the degeneracy only marginally, leav- 
ing the SN ellipses quite skinny, since the ratios /w\ in 
Figure 1 are seen to vary only weakly with z. 

2.3. Counts 

For a sample of objects volume limited out to redshift 
Zmaxj the average number per unit redshift is (W72) 



p(z) oc 



{i + z)^E{zy 



(14) 



Defining Xi — Zi, the probability distribution for the ob- 
served set of A'^ redshifts x is not a multivariate Gaussian 
as above, but a multivariate Poisson distribution, 



N 



/(x) = e-^— n5(^n), 



(15) 



where TV = (N) = /q^"*""" p{z)dz is the expected number of 
objects and g(z) = p(z)/N can be interpreted as a proba- 
bility distribution for the redshift of a typical object. Note 
that the integer N is itself random, with a Poisson distri- 
bution. Substituting equations ( p^ and ( p^ ) into (|l|) gives 



-N 



d"^ In 5 

Mm] 



g{z)dz 



1 dN dN 



(16) 



We will neglect the last term to be conservative, since it re- 
flects the information coming from the (a priori unknown) 
overall normalization. 

3. ACCURACY AND DEGENERACY 

How do these tests compare with regard to accuracy 
and degeneracy? Their degeneracy structure is illustrated 
in Figure 3, which shows contour plots of dium, df^^^/E, t, 
and D/D{0) at three redshifts. Using objects at a single 
redshift z, a test is unable to distinguish between mod- 
els lying along the same contour curve. The luminosity 
and size tests have identical degeneracy structure because 
both probe S{Kr])/K; their degeneracy curves are seen to 
rot ate anti-clockwise from a slope of 1/2 (explained below 
in § 3.1 ) at z = to negative at z — 00. The count contours 
rotate in the same sense as z increases. The isochrones ro- 
tate similarly but have a richer structure at z = because 



the age of the Universe probes E at all redshifts. They 
become vertical at high redshift where the age is inde- 
pendent of Oa. The growth factor degeneracy curves are 
seen to have a slope steeper than —1 in most of our pa- 
rameter space. This is because increasing the hyperbolic 
curvature 1 — Vim — makes fluctuation growth freeze 
out earlier, increasing D{z)/D{0), and increasing Oa at 
fixed curvature typically has the same effect. The evo- 
lution of cluster abundance places powerful constraints on 
D{z)/D{0) (Bahcall & Fan 1998). Although this test gives 
highly non-Gaussian errors e„ (the constraints are mainly 
one-sided), its degeneracy structure is still given by Fig. 3. 

The list of geometry tests that we have discussed is far 
from complete. For instance, nonlinear effects in weak 
Icnsing (Jain & Seljak 1997) and strong lensing (Falco et 
al. 1998; Bartelmann et al. 1998) are promising probes of 
r2,„ and SIa- With CMB fixing other parameters, baryonic 
features detected in future galaxy redshift surveys would 
give fairly vertical degeneracy curves, potentially measur- 
ing Vim to percent levels (Eisenstein et al. 1998). 

For all tests modeled above, the size of the error ellipses 
scales as a/\/N, whereas the shape (slope and eccentricity) 
is given by the degeneracy structure. The CMB ellipses in 
Figure 2 have been computed as in Eisenstein et al. (1998), 
marginalizing over 10 additional parameters. This CMB 
information on and J^a comes mainly from the angular 
location of acoustic features in the power spectrum, which 
depends principally on the curvature term k, i.e., on the 
combination + f^A • 
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FIG. 3 — How the degeneracy structure of different cosmological 
tests rotates with redshift. All 12 panels have same axes. 

3.1. Low redshift observations such as SDSS 

It is well known that if data is available only for z ^ 1, 
then to first order, the luminosity, angle and count tests 
probe only the parameter combination go = f7m/2 — JIa- 
In this limit, our results reduce to 



1/4 -1/2 
-1/2 1 



where Ago = 2 In lOAm/SiV^/^z 



(Ago)- 



(17) 



for luminosity tests 
with magnitude er- 



using objects at rms redshift of z. 

rors Am, Aq^ — 2(7 /N^^'^Zrms for corresponding angu- 
lar size tests on objects with fractional size errors ct, and 
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Aqo = 2{5/3Ny^'^ /zmax for number count tests volume 
limited to z^ax- This is why the corresponding z — 
panels in Figure 3 both give the same slope 1/2. 

Because of this scaling, the huge number of galaxies in 
upcoming surveys such as SDSS and 2dF may allow them 
to place competitive constraints on go, as shown in Table 
1, despite being a factor of several below SN la in redshift. 
Here we have assumed that fitting a Schechter luminos- 
ity function to N galaxies at the same redshift determines 
the parameter to within 5/-/V^/^ magnitudes, which is 
conservative based on Table 2 in Lin et al. (1996). An 
obvious obstacle to such measurements is that galaxy evo- 
lution (in luminosity, size and number density) can mimic 
a change in go- However, the brute force statistical power 
of these data sets is so large that even subsamples of 1% 
of the galaxies give interesting constraints. Studying how 
the "go" -estimates vary as the galaxies are subdivided by, 
e.g., morphology, luminosity and surface brightness there- 
fore holds the potential of providing interesting informa- 
tion about galaxy evolution and perhaps the true go . 

4. CONCLUSIONS 

In conclusion, we have derived useful expressions for the 
Fisher information matrix for a number of classical cosmo- 
logical tests and combined them with the Fisher matrix of 
the CMB. Whereas two identical data sets give only a fac- 
tor of \/2 improvement in error bars when combined, the 
gain factor was found to exceed 10 when combining SN la 
with CMB. This "cosmic complementarity" is due to the 
fortuitous fact that although cither data set alone suffers 
from a serious degeneracy problem, the directions in which 



they are insensitive (in which the ellipses in Figure 2 are 
elongated) are almost orthogonal. The complementarity 
is even more dramatic for a standard = 1, S^a = 
CDM cosmology (Tegmark et al. 1998), where a smaller 
ISW effect worsens the CMB degeneracy. 

Figure 3 shows that this complementarity is rather 
generic, with degeneracy curves in virtually all directions. 
This means that when three different tests are combined, 
there will be an important cosmic consistency check. If 
three skinny ellipses fail to overlap, at least one measure- 
ment must be wrong, whereas if they all cross at the same 
point, even hardened sceptics are likely to be impressed. 

The potential power of upcoming CMB measurements 
has led to a widespread feeling that they will completely 
dominate cosmological parameter estimation, with other 
types of experiments making only marginal contributions. 
Because of cosmic complementarity, of which the present 
paper gives a number of examples, this view is misleading: 
two data sets combined can be much more useful than 
either one alone. 
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